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Abstract 

This  report  contains  an  account  of  geostrophic  vortices 
on  a  rotating  sphere.   A  vortex  of  this  type  is  characterized 
by  a  singular  spherical  harmonic  of  degree  v  and  order  zero. 
This  function  has  a  role  analogous  to  that  of  the  Bessel 
function  which  characterizes  a  rectilinear  geostrophic  vortex 
in  a  rotating  plane. 


Ill 


1.  Introduction 

A  thin  spherical  layer  of  an  incompressible,  inviscid  fluid 
which  is  held  on  the  surface  of  a  rotating  ball  by  gravitation 
can  be  taken  for  some  purposes  as  an  approximation  to  the  Earth's 
atmosphere.   The  analysis  of  the  two  dimensional  vortical  motion 
in  such  a  layer  should  be  useful  for  the  understanding  of  certain 
observed  meteorological  phenomena. 

In  a  previous  report  a  linearized  analysis  was  presented 
under  the  assumption  that  the  outer  surface  of  the  layer  is  a 
sphere  always  concentric  to  the  ball  representing  the  Earth. 
This  report  presents  a  linear  analysis  in  which  the  outer  surface 
of  the  fluid  layer  is  allowed  to  be  free  but  the  Coriolis  force 
is  assiuied  to  be  independent  of  latitude.   This  assumption  de- 
fines a  motion  which  can  be  regarded  as  a  first  approximation  to 
the  actual  motion  even  if  the  fluid  layer  is  not  confined  to  a 
narrow  band  bounded  by  two  circles  of  latitude. 

The  investigation  covers  the  case  in  which  tangential 
accelerations  are  neglected  in  comparison  with  the  Coriolis  force. 
In  other  words,  we  study  geostrophic  vortices  on  a  sphere  as 
contrasted  with  geostrophic  vortices  on  a  plane.   The  latter  have 
been  discussed  by  several  authors  in  connection  with  a  tangential 
plane  approximation  to  the  motion  of  a  thin  layer  of  fluid 
covering  a  rotating  ball.   References  can  be  found  in  Morikawa, 

[1]. 


2.  Equations  of  Motion 

Let  p  denote  the  distance  of  a  point  from  the  center  of  a 
ball  E  of  large  radius  a  which  rotates  with  constant  angular 
velocity  to  about  a  polar  axis.   Let  ^   and  6    denote  respectively 
the  longitude  and  the  colatitude  of  a  point  on  the  rotating 
spherical  surface  S  of  E.   Let  p  =  a  and  p  =  a  +  h{^,9,t)    repre- 
sent two  surfaces  which  contain  an  incompressible^  inviscld  fluid 
which  is  gravitationally  attracted  by  E.   Suppose  h  is  small 
compared  with  a. 

The  general  problem  is  to  find  the  motion  of  the  fluid 
after  its  constant  rotatory  motion  is  disturbed  by  the  sudden 
creation,  at  a  reference  time  t  =  0,  of  concentrated  vortices 
normal  to  S. 

The  velocity  of  a  fluid  particle  relative  to  S  is  defined 
by  the  components 

u  =  (p  sin  6}  -^   =  tangential  component  toward  the  east; 

V  =       -  p  -rr  =  tangential  component  toward  the  north; 


d 

w  = 


-T&  =  radial  component, 


If  the  only  body  force  acting  is  that  due  to  the  gravitational 
potential  G  of  E,  then  the  basic  hydrodynamical  equations  which 
define  the  motion  of  the  fluid  relative  to  S  are  the  continuity 
equation 


dw  ,  ^   ,    1 


p  ^  +  2w  +  — .  ^  a 
^  dp         sm  6 


^u    ^(v  sin  0) 
5^       ^  _ 


0 


and  the  momentum  equations 

du        uv   cot   0    ,    uw    ,    „  •       a      o  „    o  1  1 

^  -  +  —  +   2cDW  sm   e  -  2cov   cos    0   =   -  ■= ^  .^    q  -^t—  > 

dT  p  p  5   p   sin  e   d^ 


c^ 


dv    ,    vw    ,    u      cot   0    ,    o  ,,   ^^o    a  11. 

-j-r-  +  —  +  +   2aou   cos    e   =  ■= -^73—  ; 

dt         p  p  5    p    d0 


2  2  1      ^Pi 

dWV  U  r^  -a  11 

4r-  -  —  -  —  -    2a)U   sm   0   =   -  ■^—  -^ —   . 
dt         p  p  ^o   ^P 

In  these   equations   the   differential   operator  with  respect   to  the 
time  means 

'^    =    ^    +  ^  ^     -  I    ^    +  w    ^ 

"HT       "5t        p    sin  0  "5^       p  "50  "Sp 

The  symbol  5   denotes  the  constant  density  of  the  fluid;  and  p 


stands  for  the  pressure  in 


^   2  2   .  2q 
5  p  CD  sxn  0 

p^  =  p  +  5^G 2 


which  can  be  referred  to  as  the  modified  pressure. 

Since  h  is  supposed  to  be  small  compared  with  the  large 
radius  a  of  E;  and  since 

w((t),  0,  a,  t )  =  0 

let  us  neglect  the  radial  velocity  and  the  radial  variation  of  u 
and  V.   Let  us  assume  that  G  =  gp  where  g  is  a  constant;  and  that 

the  centrifugal  effects  manifested  by  the  partial  derivatives  of 

2  2    2 
5  p  03   sin  0 

"    can  be  ignored.   Let  us  also  assume  that  the  motion 

is  such  that  the  nonlinear  terms  in  the  tangential  momentum 


equations  can  be  neglected;  and  that  the  radial  momentum  equa- 
tion can  be  replaced  with  the  hydrostatic  law 

pi^,G,p,t)    =   g5^(h+  a  -  p) 

which  satisfies  the  condition  that  the  pressure  is  zero  at  the 
free  surface  p  =  a+h. 

Under  the  above  assumptions  and  with  the  notation 

u{^,e,  t)  =  u{^,e,  a,  t)  , 

v{^,0,  t)  =  v{^,e,a.,t)    ; 

an  approximation  to  the  motion  is  determined  by  the  equations 


(2.1) 
(2.2) 

(2.3) 


^w 


^u 
St 


.  '\j  .   .'V. 


+ 


1     ^u   h{v   sin  e) 


a  sm 


e  m 


2(jov  COS  0  =  - 


g     ^h  . 
a  sin  6   '^  ' 


=   0  ; 


Sv      '^ 

-vr  +  2cDU  COS 


^  g  ^h 

a  "5^ 


An  integration  of  (2.1)  gives 


w(({),  0,  a+hj  t ) 


h 
a  sin  0 


S( V  sin  0 ) 


"35" 


^U 


The  kinematic  condition  at  the  free  surface  is 


w(<t),e,a+h,t)  =  ^ 


and  hence 


(2.4) 


dh  ^    h 
dT   a  sin  6 


^(v  sin  0)  _  ^u 


If  h{^6,0)   =  h     =   const.;  and  if  we  introduce 


h-  h 


r]{^,e,t)    =   -^ 


o 


assuming  that  r)  is  small  almost  everywhere,  then  a  linearization 
of  (2.4)  yields 


(2.5) 


't   a  sin 


^(v  sin  0)    ^u 
^ ^ 


while  the  momentum  equations  become 


(2.6) 


(2.7) 


■2^  -  2cDV  cos 


Sv  ,  „  ~ 

--rr  +  2a)U  cos 


^^^0 

a 

sm 

0 

a 

^5    " 

3 .  Geostrophic  Vortices 

If  we  neglect  the  variation  of  the  Coriolis  force  with 
latitude  and  take 


0)  cos  6    =    (X)    cos  0-,  =  CD-] 


the  last  equations  of  Section  2  reduce  to 


(3.1) 


^t 


a  sm 


^(v  sin  0)    ^u 
5S       "^ 


(3.2) 


(3.3) 


i  -  ^-1^ 


|I  .  2.,u 


a  sin  0  ^  ' 


gh. 


a 
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Hereafter  we  confine  ourselves  to  a  study  of  these  equations. 
As  will  be  explained  in  the  sequel,  they  lead  to  what  are  called 
geostrophic  vortices. 

The  elimination  of  t]  from  (3c2)  and  (3.3)  leads  to 


(3.4) 


St 


a  sm 


Su  sin  0  ^v 
^  h^ 


=     2C0^T1^ 


The  quantity 


1 


a  sm 


^(u  sin  0)    hv 
^  S^ 


is  the  radial  component  of  vorticity;  and  by  integration  of  (3. 4) 
we  have 


(3.5) 


a  sin  0 


h{u.   sin  0)    ^v 

— 5s         i;^ 


2a)^Ti  +i:{^,e,0) 


This  implies  that  we  can  use  (3-5)^  (3.2)  and  (3.3)  as  a  basic 
set  of  equations  Instead  of  (3.1);  (3.2)  and  (3.3). 
Equations  (3.2)  and  (3.3)  give 


(3.6) 
and 


a 

2-1 13- 

1            h^Vi 

sin  0    ^t^(j) 

(3.7) 


2~ 


ht' 


gh. 


a 
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1    OT]     d  T) 

;in  e  '^       'StdU 


If  these  are  used  to  eliminate  u  and  v  from  (3»5)^  we  find 


(3.8) 


1        ^   (sin   0)1^4-  -1^1^ 
'^^   ^    ^^  ^^         sin   0    h^^ 


2 


a 


gh 


o 


4fl>-  Tl     +    J 


fCD    Tl     +    ^ 


a^2cu^C((}),e,0) 


For  a  vortex  of  strength  \i   concentrated  at  {^^,6^)   when 
t  =  0  we  take 


'1'    !■ 


a^2ao  C  (<!>,©.  0)    a^2a)-j_     5  ( (f)  -  cf)^  )5  (6  -  0-^  ) 


"gH" 


gh. 


2    . 
a  sm 


1 


where  6  symbolizes  the  Dirac  delta  function. 

The  theory  of  Laplace  transforms  can  be  used  to  show  that 
the  steady  state  solution  of  (3.8)  corresponding  to  a  con- 
centrated vortex  is  such  that 


(3.9) 


2cjo^(({),e) 

t^oo  i^^'    '     '  gh 


o 


where  tj/  must   satisfy 


(3.10) 


.2   2 


sm 


1       4   (sin  0)   |t  +  -1-  ^  -  — p^ 
ITT  ^   ^  ^        sin^e   ^^  Sh^ 


n5((t)-  <l)i)5(e  -e^) 

sin   0- 


1 


The  function  tp   does  not  depend  on  t  and  the  associated  time 
independent  velocity  comonents  are 


:3.ii) 


(3.12) 


U=  L,_u(^,0,t)  =i|g  ; 


t^oo 


a  "5F 


V  =  L^^_v(c(),0,t) 


t^OD 


1      ^f 

a  sin  6   "5^ 


Note  that  (3.10)  is  a  consequence  of  the  equations  (3.11)^  (3.12) 
and  the  vorticity  equation 


(3.13) 


a  sxn 


^(^^^-^)-| 


^^A^oo"^  +C('l>.e.O) 
4co^^    M.5(ct)-  4^)5(0  -e^) 


gn" 


2  . 
a  sm 


1 


In  other  words,  (3.10)  is  implied  by  (3*5);  (3.2)  and  (3-3)  when 
we  ignore  inertial  forces. 

The  motion  in  a  thin  planar  layer  of  fluid  tangential  to 
the  surface  of  the  Earth  is  often  used  as  an  approximation  to  the 
actual  motion  of  the  Earth's  atmosphere.   For  such  an  approxima- 
tion the  analogues  of  (3.11);  (3.12)  and  (3. 13)  are 


(3.14) 

(3.15) 
and 

(3.16) 


u  = 


"5F 


1 
r 


SrU  ^ 


V  =  - 


4co' 


1  ^X  . 


^X   n5((j)  -  (t)^)5(r  -  r^) 

gH  r~, 

'^^  o  1 


In  terms  of  the  polar  coordinates  {r,^),    these  show  that  )(  must 
satisfy 

The  only  physically  admissible  solution  of  (3.17)  is 


(^°1Q)    X   =  -^^o 


2co 


/gh 


-  /r  +  r-j  -2rr,  cos  ( (|)  -  c|)-,  ) 


where  K  (   )  denotes  the  zeroth  order  modified  Bessel  function 

o    ' 

of  the  second  kind.  This  defines  what  is  called  a  geostrophic 
vortex.  The  motion  of  various  configurations  of  such  vortices 
has  been  studied  by  various  authors,  notably  G.  K.  Morikawa  [1], 
[2]  whose  papers  contain  detailed  explanations  and  other 
references.  In  keeping  with  what  seems  to  be  accepted  termin- 
ology we  can  say  that  (3.IO)  defines  a  geostrophic  vortex  on  a 
sphere. 

Our  object  now  is  to  investigate  the  nature  of  the  solution 
of  (3.10).   This  equation  can  be  written  as 


(3.19) 


wnere 


(3.20) 


H6(([.  -([)^)5(0  -e^) 
~  sin  6-. 


v(v+l) 


;i  2  2 
403-,  a 


The  solution  of  (3.19)  subject  to 


a   <   0    <   ^   ;         0<(f)<27r; 


and  prescrl"bed  boundary  conditions  on  the  circles  of  latitude 


e   =   a   ;         0  =  P  ; 


namely  a  Green's  function  for  a  zone,  can  be  expressed  in  the 
form 


00 

(3.21)    ^  =  ^IZ  t\^(^°^  0)+B^Q^(cos  e)]cos  m((t)-<l)-L) 

in=0 


where  F^   and  Q  are  solutions  of 

V         V 


^   (l-z2)  ^  + 
^      '  dz 


^ 


v(v+l)  - 


m 


1  -z 


f  =  0  . 


The  appropriate  series  (Jo 21)  can  be  summed  for  the  case  in  which 
the  zone  becomes  a  sphere.   However,  since  the  details  of  the 
analysis  based  on  the  use  of  (3 "21)  are  lengthy,  let  us  proceed 
in  a  different  way. 

The  general  solution  of 


^T^^(Bine)||+v(v  +  l)*. 


0 


can  be  expressed  in  the  form 


i;  =   c-,P  (-cos  e)   +c„P  (cos  e) 
^  1  V  2  V 


when  V  is  not  an  integer  as  is  the  case  when  v  is  defined  by 


v(v+l) 


ii  2  2 


or 


10 


V  =  -  i  +  iA 


where 


A  =  o 


,^2  2 


gh. 


1 


In  the  neighborhood  of  the  north  pole  P  (cos  0)  is  continuous 
and  P  (1)  =  1;  but 

V 

P   (-cos   0)    =   P  -,  (-cos   0) 

-i+iA 

00 

2          i„   -X       r                 cos    Axdx 
=  —  cosh  Att    /      

^  ^      y/2(cosh  X  -  cos   e) 


behaves    like 


0 


P  (-cos  6)   r^   .  i,n  B 

V  TT 


This  suggests  that  if  a  vortex  is  concentrated  at  {^.^,9-^);    if  d-j^ 
is  the  geodesic  distance  from  this  point  to  an  arbitrary  point 
((j),  0)  on  S;  and  if  -y-.    is  the  angle  -y-^  =  d-^/a.,    then 


fi ii P  (-cos  -v-i  ) 

4  sin  VTT  V       '^1 


4  sin  VTT  V 


-cos  0  cos 


-sin  0  sin  0-,  cos  ((j)-(|)-,) 


is  the  fundamental  solution  of  (3.19)o 

It  can  now  be  verified  by  direct  substitution  that 


(3.22) 


Th     =    -r. ti P   (-COS  7.  ) 

^   4  sm  VTT     V  '1 


does  indeed  satisfy  (3. 19).   Furthermore,  computations  with  the 
velocity  components 


11 


U  =  Ti ^ ^  P  (-COS  7-,  ) 

4a  sm  VTT  d0  V       '1 


V  = 


-n ^ 3 ^ 4r  P  (-cos  7,  )  ; 

4a  sm  0  sm  vtt  ^({)  v       '1 


show  that  in  the  neighborhood  of  {^-.,0-.),     (3.22)  defines  a 
vortical  motion;  but  the  vortex  point  itself  remains  at  rest  —  it 
possesses  no  autonomous  motion.   In  fact,  it  can  be  shown  that 


L        71 r- P  (-cos  -V-i  ) 

a->cD  4  sm  vtt  v      '1 


K 


20D 


/gh 


1   /  P    2 

-  /r  +  r-,  -  2rr,  cos  (tj)  -<})•,  ) 


Hence  we  can  say  that  (3«22)  represents  a  geos trophic  vortex  on 
a  rotating  sphere  in  the  same  way  as  we  say  that 


[I 
X  =  -  n¥ 


K 


2a) 


/gh 


—  ^r  +  r-,  -  2rr-,  cos  ( <{)  -  (|)-,  ) 


represents  a  geostrophic  vortex  on  a  rotating  plane. 

For  an  arbitrary  distribution  of  n  vortices  on  the  sphere 
the  function  ^,  which  can  be  regarded  as  a  stream  function,  is 

n 


1 

lb  =   71 : 5   U  •  P  (-cos  -v.  ) 

^   4  sin  VTT  4 — r  1  V  ^      '^x^ 
1=1 


where  -v.  =  d./a.   and  d.  is  the  geodesic  distance  on  S  from  the 
point  of  concentration  of  the  i-th  vortex,  {^.,6.),    to  an 
arbitrary  point  (({>,©).   That  is. 


cos-v.  =  cos  6   cos  0.+sin0  sin  6.  cos  (i-ct.)  . 


12 


The  associated  velocity  field  is  given  by 


U 


=  ^|l  =  4a  sin  vir   ^  III  ^^i^v  ^'^"^  1^)    i 


V  = 


a 


^   .  1^  =  ^ ^ J ^ ^^ H.P  (-cos  7.) 

sm  0   ^   4a  sm  Q    sm  vir  ^  4-rT  1  v  ^      '^i 


Since 


U  =  a  sin 


d(|) 

dt 


d0 


the  equations  of  motion  of  the  vortex  at  (<[),  ,6,  )  are 


and 


d^ 


I,  2  . 
4a  sm 


sm  VTT 


n 


4>=4, 


dS^ 

dt" 


4a  sm 


sm  vir 


n 


^IZ^^iPv^-^^^  ^i^ 


0=0 


¥e  are  now  in  a  position  to  develop  a  general  theory  for 
geostrophic  vortices  on  a  sphere  analogous  to  the  basic  Helmholtz- 
Kirchoff  theory  for  rectilinear  vortices  normal  to  a  plane. 
Instead  of  doing  this  here  it  seems  more  useful  to  examine  a  few 
special  cases  which  may  be  approximately  applicable  to  certain 
physical  situations. 


13 


4.  Geostrophic  Vortices  in  the  Northern  Hemisphere.   Either 
Velocity  Component  Zero  at  the  Equator 

Let  us  turn  now  to  the  case  of  a  concentrated  geostrophic 
vortex  in  the  northern  hemisphere  subject  to  the  condition  that 
the  normal  velocity  of  the  fluid  at  the  equator  is  zero.   For 
this  case  the  stream  function  tjj   must  satisfy 

sm  0  dtp 

sin  U-. 
for 

0    <   0,e^    <  5  ;         0   <   (j),  (j)^    <   27r    . 
The   boundary   condition   is 
(^.2)  V((j),l)    =i^^(4,,|)    =    0 


or 

200 


^  ^((f.,|)    =   L^^^ii((t),0,t)    =   const, 


We  have  supposed  that  ri{^,—,0)    =   0;    and  we  now  suppose  that 
T]{^,^,t)   =   0  for  all  t.   The  boundary  condition  (■^.2)  then  be- 
comes 


(^.3)  i^i^,^)   =  0  . 

The  solution  of  (4.1)  which  satisfies  (4.3)  can  be  formed 
from  the  fundamental  solution 


14 


n y: P  (-cos  -v,  ) 

4  sin  VTT   V       '^l 


is  we  assume  that  the  idea  of  reflection  across  the  equator  will 
lead  toward  the  solution.   In  fact,  it  is  easy  to  verify  that 


/^P[-cos0cos0,  -sinSsinS,  cos  (c|)-(j)-,)]; 

(4.4)     tp  =  J, — ^ —  J  y 

^     ^     ^    4  SXn  VTT  \  ( 

-P  [cos0cose,  -sin0sin0,  cos  ((f)  -(f)-,)] 


is  the  desired  solution.   With  respect  to  the  whole  sphere,  ^  of 
(4.4)  has  concentrated  vortices  at  ((j)-,,0-,)  and  ((j)  ,Tr-0-,).   It 
should  be  noted  that  if  we  seek  a  stream  function  for  a  geo- 
strophic  vortex  on  a  rotating  plane  such  that  the  normal  velocity 
of  the  fluid  is  zero  along  a  circle  whose  center  is  at  the  center 
of  rotation,  then  a  result  analogous  to  (4.4)  does  not  exist. 

The  velocity  components  of  the  motion  defined  by  (4.4)  are 


U  =  - 


1  ^ 
a  "SS" 


% 


V. 


a  sm  VTT 


-cos  Q  cos  Q- 


-sin  0  sine,  cos  (({)-(j),) 


sin  Q  cos  e 


1 


'+P" 


cos  0  cos 


-sin  e  sin  6-,  cos  ((|)-  (|)  ) 


-cos  Q  sine,  cos  ((|)-<|5-,  ) 

sin  e  cos  e, 
+COS  e  sin  e  cos  (({)-  ({)-,  ) 


and 
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V  = 


a  sin  6    ^ 


^ 


^i 


-cos  Q   cos 


a  sm  vfT 


-P' 


-sm 


1 


-sin  0  sin  0-,  cos  {^  -  (f)-,  ) 


cos  0  cos 


sin  0-.  cos  ((j)  -  (f)-,  ) 


sin  0-,  sin  ((|)  -  (i"-]  ) 


.  sin  0-,  sin  {^  -  \-.) 


The  equations  for  the  motion  of  the  vortex  concentrated  at 

[^   ,e    )    come  from  evaluating  at  {^-.,9    )    the  velocity  field  due 

to  the  vortex  at  {^-^,tt-Q^)  .      They  are 


^1  = 


dT 


cos  0 


2         V 
2a  sin  vtt 


i —  p' (cos  20.  )  ; 
V       1  ' 


=  0  . 


From  these  we  see  that  the  vortex  remains  on  its  initial  circle 
of  latitude  which  it  traverses  with  constant  angular  velocity. 
The  function  ip   corresponding  to  a  concentrated  vortex  in 
the  northern  hemisphere  with  the  boundary  condition 

U(^,^)  =  0 


can  also  be  found  by  appealing  to  a  reflection  process.   It  is 


^   = 


r 


kL 


sin    VTT 


F   [-cos    0    cos    0-,  -  sin   0    sin   0-,    cos    (tf)-  \>j)] 
+P   [cos   0   cos   0      -    sin  0    sin  0-,    cos    (<{)  -  (l)-j_)] 
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It  should  be  remarked,  however,  that  a  reflection  technique  for 
the  determination  of  solutions  of  (5.19)  for  a  general  zone  is 
unknown. 

If  the  northern  component  of  the  velocity  of  the  fluid  is 
zero  at  the  equator;  and  if  we  have  a  polar  vortex  of  strength  \i 
with  three  others  of  equal  strength  ^l  which  at  t  =  0  are  situated 
at 


4 


(^i,0i)  ;  {^  +  ^,0^)   ;  {^^  +  tlL,e^) 


in  the  northern  hemisphere,  then  the  stream  function  is 


i^   =  Tf 


1 


sm  vir 


\i    [P   (-cos   e)  -P   (cos   6)1 


+[l 


n2Tr. 


'2_   r  V   [-cos  e  cos  0-^  -  sin  6  sin  6     cos  {^-^-.  -  ^^] 


n=0 


n2Tn 


-P  [cos  6  cos  e       -  sin  Q  sin  0.,  cos  ((|)-(f)-,  - --|^]  7 


A  calculation  shows  that  the  polar  vortex  remains  stationary 
while  the  paths  of -the  other  three  are  determined  by 


0v  =  0 


r      n^[P^(-cos  6^)   +P^(cos  e^)]  I 


(^•5)      k 


4a  sm  VTT 


+M. 


,'^3 


3  cos  0iP^(f  sin^e^-  1) 

2  cos  0-,p'  (cos  20-,  ) 
1  V       1^ 

+  cos  0iP^(l  -^  sin^e.^) 


O 
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In  other  words,  the  position  of  the  vortices  subsequent  to  t  =  0 
are  given  by 

(<t,^+Ot,0-^)  ;    (^^  +  ^+Ot,e^)  ;    ((f.^+^+Ot,0^)  . 

Equation  (4.5)  shows  that  for  a  certain  strength  \i      the  angular 
velocity  O  becomes  zero  and  all  the  vortices  remain  stationary. 
The  stability  of  this  equilibrium  configuration  can  be  investiga- 
ted by  using  the  methods  explained  in  the  paper  by  Morikawa  [1]; 
and  the  paper  by  Morikawa  and  Swenson  [2]. 
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